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We report on the valley blockade and the mnltielectron Kondo effect generated by an impurity 
atom in a silicon nano field effect device. According to the spin-valley nature of tunnelling processes, 
and consistently with those allowed by the valley blockade regime, the manifestation of Kondo effect 
obeys to the periodicity 4 of the electron filling sequence typical of silicon emerging at occupation 
N — 1, 2, 3. The spin-valley Kondo effect emerges under different kinds of screening depending on the 
electron filling. By exploiting the valley blockade regime, valley index conservation in the Kondo 
SU(4) is deduced without the employment of an external magnetic field. Microwave irradiation 
suppresses the Kondo effect at occupancies up to three electrons. 


I. INTRODUCTION 


The control of individual electrons and impurity atoms 
in silicon nano Field Effect Transistors (FETs) has deep 
impact on the field of valleytronicJ^I, which consists of 
exploiting the orbital part of the electron wave function 
as additional degree of freedom with respect to more 
conventional charge and spin states. Valley-related 
effects have been lately reported, from the vallw Riling 
sequence in silicon quantum dots (Si-QDs)2ffil and 
lifetime-enhanced transport as indirect observation of 
spin and valley blockad^, to a tunable valley Kondo 
effect in a As atom at single electron filling (N = ijSl. 
Valley-based qubit^ and pure valley blockad^ have 
been proposed in analogy to spin qubits and Pauli spin 
blockade respectively. The direct experimental obser¬ 
vation of the valley blockade, namely the suppression 
of transport determined by the orthogonality of orbital 
degrees in a reservoir and at an impurity site, was till 
now unaddressed. Valley quantum numbers may also 
influence higher-order tunneling processes involving two 
or more particles, as in the case of the Kondo effect. The 
simplest Kondo theory for III-V semiconductor quantum 
dots accounts for the emergence of Kondo effect only 
with uncoupled spins in the doiP^ii^. Generalized to the 
case of Si-QDs and group-V donors, the Kondo physics 
predicts that both spin and valley indices of a confined 
electron can be screened by a Fermi sea through virtual 
spin-valley flip processes, allowing for highly symmetric 
SU(4) Kondo stateJi^. 

We already demonstrated single/few donors based silicon 
transistors operabilitjli^lR^ and microwave irradiation 
effects in single atom based devices for spectroscopy of 
excited states at zero biaJi^ as well as in single interface 
defect^am. 

Here we investigate a system of a P atom strongly 
coupled to the leads of a silicon ultra-scaled tri-gate 
FET (Fig. Peculiarities of valley dependent 

quantum transport arise both in sequential tunneling 
processes and in Kondo-like cotunneling events. The 


amount of valley mixing in contact-impurity tunnelings 
gives rise to different kinds of Kondo screening from 
iV = 1 to V = 3 at the relatively high temperature 
of 4.2 K. In agreement with theory, standard Coulomb 
blockade is observed at V = 4 where Kondo effect is 
forbidden. The fourfold degeneracy is usually revealed 
by splitting the orbital scr eening from the spin screen¬ 
ing with a magnetic fielcpEnMl, Differently, in our 
experiment the valley blockade of first-order transport 
{N = 0 —>■ 1 —>■ 0) influences the Kondo processes within 
the N = 1 Coulomb diamond: it selects donor-lead tran¬ 
sitions between identical valley states, which beautifully 
confirms the selection rule at the base of the SU(4) 
symmetry. The Kondo perturbed regime is obtained at 
the base temperature at filling up to 3 electrons and it 
is suppressed by means of microwave irradiation. 

The paper is organized as follows. In Sec. II the 
device and the experimental methods are described. 
In Sec. Ill the experimental results are presented and 
discussed. The Appendix connects the spin-valley states 
mentioned in Sec. Ill with the established formalism of 
Ref. [211 In Supplemental Material details about device 
electrical characterization and data analysis procedure 
are reported. 


II. DEVICE AND METHODS 

The tri-gate FET is fabricated from fully depleted 
silicon-on-insulator (FDSOI) technology. The silicon de¬ 
vice layer is phosphorous implanted (10^® cm“®) and then 
etched to define the nanowire channel (Fig. The gate 
is then formed on the three sides of the narrow channel 
with a 5 nm thermal silicon dioxide of isolation. Silicon 
nitride (Si 3 N 4 ) spacers around the gate protect the ac¬ 
tive region from Arsenic implantation of highly doped 
(metal-like) source-drain electrodes. Nominal channel 
width of the tested device is 50 nm, gate length is 20 
nm and channel is 8 nm thick. A dc voltage Vg sets the 
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FIG. 1. (Color online) (a) Scanning Electron Micrograph 
image of a tri-gate FET identical to measured device. Be¬ 
low the sketch of the active region with the two conductive 
paths. Brown denotes the gate contact whereas Si 3 N 4 spacers 
around the gate are yellow. Heavily-doped regions of reser¬ 
voirs are highlighted in blue, (b) From the stability diagram 
two Coulomb blockade patterns appear, each with its own 
lever-arm factor a. Red-orange (high current) diamonds are 
attributed to the donor Dqd strongly coupled with leads. Blue 
(low current) diamonds are relative to a weakly coupled quan¬ 
tum dot. 


gate potential to tune the number of confined electrons, a 
bias Vsd controls the difference between the Fermi levels 
of the electron reservoirs at the contacts. The differen¬ 
tial conductance dl/dVsd is measured by using a stan¬ 
dard lock-in excitation of 40 fJ.Vrms at 116 Hz applied to 
the source electrode. The microwave line consists of a 
3.5-mm-diameter beryllium in stainless-steel coaxial line 
(UT-141) and it is ended by an unmatched dipole an¬ 
tenna. 

As in similar dopecP^ and undopecP^ devices with tri¬ 
gate geometry, the onset of quantum transport takes 
place via two parallel paths at the corners of the 
nanowire. Two sets of peaks differing by orders of magni¬ 
tude in conductivity provide distinct capacitive couplings 
with the gate (Figs. and &)• 

The highly conductive path Dq^ is attributed to a donor 
close to the edge of the silicon nanowire with electronic 
states hybridized with the interface states of the cor¬ 
ner. For the exhaustive discussion see Supplemental 
MateriaP^, in agreement with va lues of donors embedded 
in etched-silicon channeld^^^SIlIl j^;g |;.y^ro (subthresh¬ 
old) peaks Hqjj and D~^ are associated to transport via 
neutral and negative charge states of the phosphorous 
atom respectively. At higher gate fields the potential 
barriers at the ends of the tiny channel originate quan¬ 
tum dot-like states visible as above-threshold resonances 
and Coulomb diamonds with N > 3 electrons confined. 
Such donor based quantum dot system is referred as the 
donor hereafter. 


III. RESULTS AND DISCUSSION 
A. Valley blockade 


The valley blockade arising at single electron filling is 
observed by the asymmetry of the stability diagram of 
^ T>qd transitions in Figure [^. The lack of con¬ 
duction in the section of the map involving the ground 
state only at negative bias voltages is the hallmark of 
valley blockade regime predicted in Ref. [8] as reported in 
Figure It is determined by perpendicularity of the 
orbital states at the Fermi energies of the contacts and 
the ground state of the donor, which is labeled with the 
odd parity index o. The first excited state of the donor, 
labeled with the even e parity index, defines in addition 
the inelastic cotunneling region delimited by the hori¬ 
zontal dashed lines. The valley splitting A is 6.2 me\^. 
The blockade is lifted at Vsd = —6.2 mV when the donor 
excited state is resonant with the e states at the Fermi 
level in the left lead. 

The origin of the valley blockade regime is ascribed to 
the random graining of the As atoms of the contacts in 
proximity of the active region, which originates singulari¬ 
ties in the reservoir density of stateJ^. As the valley par¬ 
ity of the wave functions depends on low dimensionality 
of randomly diffused dopant distributiorPS, the electron 
states at the Fermi levels of source and drain result un¬ 
intentionally of opposite valley. 

Such very special condition is confirmed by the stability 
diagram of the weakly coupled quantum dot in Figure 
j^. The same asymmetry as the ^ 1?°^ resonance 
proves from an independent measurement that the valley 
blockade occurs because of the opposite valley polariza¬ 
tion of the contacts. 

The effectiveness of the valley blockade over the compet¬ 
ing Coulomb blockade is quantified through a simple rate 
equation model: tunneling events between identical val¬ 
leys turn out to be more than 10 times faster than those 
between states of different valley. The current across a 
potential barrier, for example the left one, in the sequen¬ 
tial regime is Aeft = e(ri„rout)/(rin TFout), where Fin 
and Tout denote the rates in and out of the donoil^. Any 
temperature dependence is neglected since we are inter¬ 
ested in a ratio between tunneling rates and not in their 
absolute values. When n levels of the donor enter in the 
bias window, Fin has to be replaced with 
ferring to the band diagrams (b) and (c) of Figurej^, we 
define FA as the rate of tunneling events between levels 
with opposite parity (i.e. from the left lead to the donor 
ground state). Analogously, FPq corresponds to tunneling 
processes preserving the valley index (i.e. from the left 
lead to the donor excited state). In particular, along the 



spectively. Since the valley blockade selects only tunnel¬ 
ing events entering the donor, relaxations are admitted 
in the traversing time, giving to a unique outcoming rate 
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FIG. 2. (Color online) (a) Stability map of the 73°^ 

resonance and first Coulomb diamond {N = 1). Dashed lines 
indicate the onset of cotunneling via the excited state. For a 
better visualization the derivative of the differential conduc¬ 
tance is computed, (b) Expected stability map in presence of 
inelastic cotunneling and valley blockade and, on the right, 
energy diagrams for sequential transitions. Thin dashed lines 
represent empty states available in the reservoirs, (c) Relative 
comparison of tunneling rates between states with identical 
and opposite valley composition. Inset: current stability dia¬ 
gram from which the value /a and /<> used in Equation 0 are 
obtained, (d) Transport data before the onset of conduction 
through the donor. The sensitivity is tuned to magnify the 
current contribution of QD. The first peak shows the asym¬ 
metry due to the opposite valley polarization of the source 
and drain electrodes and the consequent valley blockade. GS 
and ES denote the ground and the excited state respectively. 


Tout- 

At the stationary state /left = /right = /■ Because of the 
valley blockade, we consider PA «C By further as¬ 
suming similar couplings between the donor and the two 
reservoirs (PPj ^ Pout)) we obtain: 

~ 1 __ ( 1 ) 

/O l + 2Pi^/PP„ 

/a and /<) are evaluated from the data plotted in the 
inset of Figure]^, so that the ratio PA/PPj is extracted 
for different values of gate voltage. 



T(K) 


FIG. 3. (Golor online) Gonductance at different temperatures 
spanning the first spin-valley shell, for reference see Fig. 0^. 
The increase of conductance at low temperatures is the finger¬ 
print of emerging Kondo physics. Three Kondo temperatures 
Tki = (8.2±1.3) K, Tk2 = (3.0±0.4) K and Tks = (3±1) K 
are obtained at A = 1, 2, 3 electron filling respectively. In the 
Kondo perturbed regime the fitting funct ion adop ted is suit¬ 
able for both SU(2) and SU(4) frameworkfpil^^ For A = 4 
the Kondo trend is not observed as the shell is complete. 


B. Kondo effect 


We now focus on the experimental investigation of 
the strong coupling regime to the leads in the first sil¬ 
icon spin-valley shell (i.e. from A = 1 to A = 4) of 
the donor-dot system. In III-V semiconductor quantum 
dots the presence of the Kondo effect stems from the 
paired-unpaired spin filling sequence at odd occupation 
number^SEIl differently, here the Kondo periodicity is al¬ 
tered by the presence of valley degrees of freedom. A sili¬ 
con system with an orbital spacing larger than the valley 
splitting, such as an isolated donor, is predicted to have 
a Kondo periodicity of foui^^, as in carbon nanotube^^. 
Figure 0 provides the first experimental signature of a 
spin-valley Kondo effect reflecting the fourfold degener¬ 
acy of the first orbital shell. The hallmark of the Kondo 
perturbed transport lies in the increase of conductance 
by lowering the temperature down to 4.2 K at A = 1, 
2 and 3. Such fingerprint is not observed at A = 4 as 
no unpaired spin-valley degrees can be screened by the 
Fermi sea of the leads. In Figure]^ the Kondo resonance 
at Vsd = mV is marked with a green rhombus. The 
non-zero conductance background from |I4j;| > 5 mV is 
attributed to inelastic cotunneling. To discern the con¬ 
tribution to conductance due to cotunneling from that 
Kondo-related, the opposite temperature dependence of 
the two second-order transport mechanisms is exploited. 
By heating the sample the Kondo resonance progressively 
damps till vanishing at 24 K, while the cotunneling back¬ 
ground is enhancecpS. 

The neutral charge configuration of the donor //^^ is 
shown in detail in the right part of Figure]^. The selec¬ 
tion rule on valley parity underlying the valley blockade 
regime is reflected in second-order Kondo transport: only 
those processes for which the parity index is conserved 
clearly stand out from the cotunneling background. Fig- 
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FIG. 4. (Color online) Conductance as a function of bias in 
the = 1 and N — 2 configurations. On the right relative 
Kondo processes, (a) If the valley parity is conserved the 
Kondo resonance is observed {Vsd = —7) mV, whereas at 
Vsd = 0, +7 mV the Kondo events are prevented by the valley 
blockade, (b) Schematics describing spin-SU(2) Kondo effect 
in presence of valley mixing. 


ure|^ reports a section of the Coulomb diamond. Kondo 
events with and without valley flips of the electron bound 
to the don or (exp ected at Vsd = +7 mV and 0 mV 
respectively^SEMi]^ marked by black arrows) do not ap¬ 
pear. 

The general Anderson Hamiltonian including valley de¬ 
grees of freedom reads a^ 


H = 


km<7 




-I- 


mcr 






mm' kmcr 

+ E^-v( ^\ma^rniT + d\j^^Ckma) 

km<7 


dlna-^kma) 


( 2 ) 


where , c create and annihilate non-interacting fermions 
in the leads, Sma is the single-particle energy level of 
the donor localized state with valley m and spin cr, d'' 
and d are the creation and annihilation operators of 
this state, Umm' the intra {m = m') or inter (m ^ m!) 
valley Coulomb repulsion and n = d^d is the number 
operator; c^d and d)c describe tunneling events between 
a spin-valley state in the reservoir and a spin-valley level 
of the donor and viceversa. The channels involving the 
coupling of the donor with the reservoirs are weighted 


by Vo for tunneling amplitude of processes preserving 
the valley index and Vx for exchanges between different 
valleys (m and fh indicate opposite parity indices).!^ 
The lack of a zero-bias Kondo peak reveals an intervalley 
coupling Vx smaller than the intravalley term Vq . Kondo 
processes at Vsd = — 7 mV preserve the electron valley 
parity: the Hamiltonian in Equation (H) reduces to the 
high symmetric form with just the Vq intravally coupling 
term. It corresponds to a SU(4)-symmetric Kondo model 
arising from a twofold spin degeneracy combined with 
a twofold valley degeneracy; as the cotunneling onset 
marks a non-zero valley splitting weakly dependent on 
14 (Fig. ^), the condition of nearly-degenerate valleys 
is here achieved by setting a bias equal or larger than A. 
Spin-valley Kondo usually associates with spin -valley 
blockade, as reported for carbon nanotube d^^FO l where 
reservoirs and dot are formed within the same tube. 
Here, as the wave functions in the leads are built from 
electron states of As with same orbital symmetry as 
the electron state at the P in the channel, the exper¬ 
iment reveals that in the device under investigation 
such maintenance of symmetry is achieved thanks to 
the specific distribution of As atoms in the reservoir 
areas. In HI-V semiconductors SU(4) symmetric Kondo 
emerges either from spin-charge entanglement in two 
electrostatically coupled quantum dots^S or by inducing 
multiorbital ground states in vertical dot^. In our case, 
two important clues point towards the SU(4) symmetry: 
firstly, at operating temperatures of the order of the 
Kondo temperature T^i the SU(4) physics dominates 
over the SU(2); secondly, the energy difference between 
the Kondo peak and the excited state 0.8 meV) is 
comparable to Tk-^- The combined emergence of valley 
blockade and SU(4) Kondo effect leads to a unified 
picture of quantum transport based on valley parity 
conservation during tunneling. 

Differently from N = 1, the Kondo effect at A = 2 
emerges as a zero bias resonance. The participation of 
the excited state e would imply a Kondo temperature 
Tx 2 ~ where ks is the Boltzmann constant, 

yielding Tx 2 ~ 70 K, a value that can be easly excluded 
thanks to the experimental data reported in Figure 
Therefore the single-particle level picture has to 
be replaced with two-electron states obtained from 
Slater determinants. Following the formalism previously 
adopted, the orbital part of such determinants is a 
combination of e and o states, and the spin component 
is either a singlet \S) or a triplet |T). Their total 
antisymmetric tensor products provides the two-particle 
states: |S')|e, e), |5')|o, o), |T)(|o, e) — |e, o))/E and 
|S')(|o,e) -I- |e,o))/v^. The electron-electron interaction 
further mixes the orbital components associated to 
the three spin singlets. As a result, both singlets and 
triplets have as orbital counterpart a linear combination 
of different valley parities (see Appendix). The Kondo 
transport takes origin when an electron in the reservoir 
with a well-defined valley parity tunnels through a 
barrier by changing its valley index to form the com- 
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binations of mixed valley eigenstates above mentioned. 
For this reason these Kondo exchanges do not suffer any 
effect of the valley blockade previously discussed: the 
Vx mixing channel governs the Kondo transitions within 
the N = 2 diamond. The exchange integral between 
wave functions with different combinations of valleys 
is n egligibl^^ in the limit of non-interacting spins \S) 
and |T) are degenerate in energy. Hence, each level 
admits spin flip Kondo processes at zero bias (sketch 
in Figure &)■ The ground and the first two-electron 
excited level are coupled with the conduction electrons 
of the leads through two orthogonal scattering channels, 
leading to a two-level SU(2) Kondo effect. That one 
with higher Kondo temperature is activated here. The 
vanishing of the exchange interaction hampers a univocal 
labeling of our Kondo at = 2 as under, double-fully 
or overscreened type. 

Such an anomalous Kondo effect at even occupancy 
is inherently related to the non-zero momentum of 
silicon conduction band edges. In III-V direct-gap 
semiconductor devices integer total spin states may 
lead to the appearance of Kondo peaks at even electron 
fillings. In Ref. [Ml a magnetic field brings singlet and 
triplet into degeneracy, giving rise to singlet-triplet 
Kondo transitions at zero-bias. Alternatively, by means 
of an electric field the orbital spacing of two adjacent 
levels can be tunecP^ to promote the spin triplet state 
as the ground state. In the system under investigation 
such condition is experimentally ruled out, as by varying 
Vg the Kondo resonance remains pinned at zero bias. 
Within the N = 3 diamond the contribution to 

zero-bias conductance due to the Kondo effect (third 
panel of Fig. and Fig. is attributed to the 

unpaired spin of the third electron fully screened 
the many-body spin of the electrons in the reservoirgSHli] 


C. Microwave suppression of Kondo effect 

Finally the attention is turned to the perturbation of 
the system by means of a tunable microwave irradia¬ 
tion. In this section the experimental microwave sup¬ 
pression of the Kondo effect jpreviously observed in di¬ 
rect bandgap semiconductorPSHini^ is experimentally ad¬ 
dressed in a multi-valley semiconductor like Si. 

As the setup is equipped with an unmatched coaxial 
cable employed as an antenna, the random microwave 
field distribution induces oscillations both to bias and to 
donor chemical potential. An empirical conversion be¬ 
tween nominal power P and effective amplitude of the 
perturbation 14; on the donor sit^^ is established for 
each electron occupancy (A^ = 1,2, 3) (details reported in 
Supplemental Material^. However, the microwave cou¬ 
pling with the quantum system, as well as the line trans¬ 
mission efficiency, varies also at different frequencies. 
For this reason, we first developed a robust background 
subtraction procedure to isolate the Kondo resonances 


from thermal-activated or photon-assisted backgrouncP^. 
Next, the Kondo resonance suppression is investigated 
by varying the frequency. The zero bias resonance at the 
lowest filling at which it is observed, namely N = 2 at 
Vg = 34: mV, is studied under microwave irradiation at 
different frequencies. Above 30 GHz the photon energy 
approaches kBTx 2 and the microwave-induced suppres¬ 
sion curves collapse over a single trencP^, revealing that 
in this regime the microwave field perturbs the coherent 
spin flip events underlying the Kondo transport. Because 
of the high working temperature an appreciable suppres¬ 
sion is achieved when many photons cooperate to bring 
the Kondo state out of equilibrium, eVL > hv. We under¬ 
line that such suppressive trend stands out independently 
from the choice of the background fitting functiorP^. 

The systematic comparison among the Kondo reso¬ 
nances at N = 1,2,3 by varying the power of a 30 
GHz excitation is carried out. In Figure the ratio 
{dl/dVsd)peak/[dl/dVsd)diirk Is plotted as a function of 
the parameter eV^/ksTK- Here ’’dark” corresponds to 
the total absence of any external microwave signal. Such 
renormalization allows a direct comparison of the three 
cases N = 1,2,3 taking into account for the different 
Kondo temperature. Microwave irradiation affects the 
Kondo resonance for all the three occupancies. Such 
suppression can be ascribed to the quenching effect of 
the microwave photons rather than to heating, that is 
excluded by the estimation of the electronic temperature 
by means of Coulomb blockade thermometry. Under the 
assumption of temperature-broadened peak, i.e a cosh“^ 
lineshape, we register an increase of ~ 0.8 K at the maxi¬ 
mum power (-5 dBm nominally) with respect to the base 
temperature in dark conditions. The weak thermal sup¬ 
pression expected in this range of temperatures cannot 
account for the complete destruction of the Kondo reso¬ 
nance, which has to be attributed to the microwave field. 
At amplitudes of energy oscillations smaller or compa¬ 
rable to the Kondo binding energies, such behavior is 
ascribable to adiabatic variations of the b ias co mbined 
with weak modulations of the gate potentiaP®^, though 
their separate contribution is difficult to evaluate. By 
increasing the power of irradiation the rate of photon- 
induced events increases as well, progressively inhibiting 
the Kondo processes. 

The response to microwave signal differs among the cases 
N = 1,2,3, as shown in Figure The slope of the 
suppression curve for V = 1 is smaller than that for 
N = 2, 3, which by contrast are comparable. An exten¬ 
sive theorjEHm has been developed for a single spin 1 /2 
in the Kondo regime (i.e. T Tk) perturbed by photons 
with hv > ksTx- Though a well-established theoretical 
analysis is still lacking for the mentioned working condi¬ 
tions (multi-valley semiconductor, hv ^ k^Tx) we can 
infer valuable information from our experimental data. 
The different suppressive rates of Figure [^suggest a con¬ 
nection between the Kondo symmetries of the three oc¬ 
cupancies. It is notable that the occupancies with similar 
damping rates have comparable Kondo temperatures and 
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FIG. 5. (Color online) Normalized Kondo peak amplitudes for 
N = 1,2,3 as functions of eVuilksTK for the three different 
occupancies. 



FIG. 6. (Color online) (a) Differential conductance of the Dqd 
resonance without microwave irradiation, (b) Same stability 
diagram measured at -17.5 dBm, -10 dBm, -7.5 dBm, -5 dBm, 
-2.5 dBm and 0 dBm nominal powers with a frequency of 
30 GHz. The asymmetry due to the valley blockade is not 
influenced by the microwave irradiation. 


the same underlying symmetry. Such suppressive trends 
are qualitatively in agreem ent w ith previous observations 
in SU(2) spin Kondo effeclP^*^. 

The reasons of the different rate for = 1 are not un¬ 
ambiguously addressed. On the one hand, the microwave 
suppression could be less efficient on SU(4) than on SU(2) 
symmetry due to different couplings of the microwave 
field to spin and to valley degrees. On the other hand, 
as Tr -1 > stating the different symmetry at 

= 1, the energy carried by 30 GHz photons may result 
to be less effective in perturbing the Kondo effect because 
of the higher binding energy ksTKi- 
Overall, our results show the effectiveness of microwave 
suppression over the first spin-valley shell, suggesting a 
possible extension of the theory of Refs. El El beyond 
the context of standard spin 1/2 Kondo effect. 

Finally, we investigate the effect of microwave irradiation 
on the valley blockade regime. The valley selection rules 
are not perturbed by the 30 GHz microwave irradiation. 
In Figure no significant variations with respect to the 
dark condition appear in the asymmetric 
peak by increasing the microwave power. At the Fermi 
levels of the reservoirs the valley parity electronic states 
are not altered and the valley blockade phenomenology 
is fully preservecP. Such experimental result suggests 
that the microwave suppression of the Kondo resonance is 
mainly due to spin coherence breaking, while valley fluc¬ 
tuations remain unaffected. The Kondo peak at A^ = 1 is 
effectively suppressed at VL = 2 m\^. Non-zero peaks 
of Figure can be therefore ascribed to competing deco- 
herent mechanisms like thermal fluctuations or inelastic 
cotunneling. 


CONCLUSIONS 

Valley blockade and multielectron Kondo related trans¬ 
port are observed in a single P donor silicon transistor. 
The Kondo perturbed regime arises at partial filling of 
the first spin-valley shell of silicon, V = 1, 2,3. At V = 1 
occupancy the pure valley blockade occuring between 
donor and leads selects Kondo spin-valley processes at 
non-zero bias consistent with a SU(4) symmetry. Mixed 
even-odd valley states originate spin SU(2) Kondo reso¬ 
nances at zero bias when 2 and 3 electrons are confined. 
The Kondo effect is not observed when the shell is com¬ 
plete at Af = 4, as predicted. Finally, microwave driven 
oscillations suppress all the three Kondo resonances while 
leaving the valley blockade phenomenology unaffected. 
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APPENDIX 

In this Appendix the spin-valley states are explicitly 
expressed. In Table we connect our formalism to the 
effective-mass theory of Hada and EtcP^ by writing the 
states in both bases. For the two-particle states the su¬ 
perscript 1 labels the first electron wave function and 2 
the second electron wave function. 
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By taking into account the Coulomb interaction 
4 Ti-eeoki-r 2 | t>etween the two electrons as a perturbation, 
the eigenfunctions can be written as linear combination 
of Slater determinants. The Coulomb interaction is spin 
independent, so mixed singlet-triplet states annihilate, 
and just combinations between the three singlets, and 


the triplet apart, remain. An example of the orbital 
part of a spin singlet is [—i sin20(|e, e) + |o, o)) + 

cos20(|o, e) -I- |e,o))], which shows how both even and 
odd valley parities contribute. 



Spin 

Orbital 



Even-odd formalism 

Ref. 122] formalism 

Single 

particle 

1 i>, 1 t) 

|e) 

|o) 




|ee) 



|5) 

|oo) 


Two particles 


^(|eo) -b |oe)) 


|F+) 




IF-) 

^(|eo) - |oe)) 



|Fo) 




TABLE I. Conversion table between the even-odd formalism adopted in the main text and the formalism of Ref. [22] 


SUPPLEMENTAL MATERIAL 
Estimated energies of the P atom 

In Figure lb of the main text two series of conductance peaks differing by an order of magnitude appear. Red-orange 
(i.e. high-current) diamonds give a unique lever-arm factor a, thus a single coupling with the gate, observed up to 
N = 5. The blue diamonds of Figure lb have a different lever-arm factor. The pair of interlaced diamond patterns 
with different lever-arm factors suggests a circuital scheme with two conductive paths in parallel (Dqd and QD) with 
a negligible inter-conductanc^^. The first is associated to a phosphorus atom, labeled as Dqd, while the second to a 
disorder-assisted quantum dot (QD), like previously observed in similar samples. Dqd localized states arise from the 
hybridization of Coulomb potential of a donor, with the gate-induced potential well at the Si/Si02 interfac^^. 

According to the doping concentration, we expect 6-12 donors in the volume of the channel. By lowering temperature 
thermal transport is progressively quenched so that the threshold voltage Vth can be easly extracted from the 
linear part of the G{Vg) characteristics, see Figure At 300 K the subthreshold current is ascribed to thermally 
broadened tunnelings through the dopants of the active regiorP^ (Fig- At cryogenic temperatures only those 

donors sufficiently coupled with both source and drain may be observed from quantum transport. Similarly to previous 
reports on single atom transistors with etched channelJi^ElIlIl^ the first addition energy of Dqd lies in the range 15 — 40 
meV (here U = 26.4 meV). The neutral charge state is 44.7 meV below the conduction band edge marked by 
the threshold voltage. This ionization ener gy is very close to the 45 meV of P atoms in bulk silicoiP^, indicating 
the proximity of the atom to the gate oxide^. Such a vicinity enhances the confinement in the gate-field direction, 
delocalizing the electron wavefunction along the tunneling direction. The consequent strong tunnel coupling with 

















FIG. 7. (a) The threshold voltage of the Dqd conductive path is extracted from the linear part of the GiVg) characteristics at 
relatively high temperatures. The contribution to G of the QD island is here neglected being one order of magnitude weaker, 
(b) The presence of donors in the channel is confirmed by the subthreshold current at room temperature. 


a) 


b) 


Background: 2-Boltzmann 




Vsd (mV) 


Vsd (mV) 


FIG. 8. (a) Raw data of the Kondo resonance for = 2 at different microwave power at 30 GHz frequency. Dashed arrows 
identify the background enhancement due to photon assisted tunneling, whereas the solid arrow highlights the suppressive 
trend of the Kondo resonance, (b) Same data after background subtraction. 


reservoirs allows for cotunneling and Kondo processes. 

The conclusive argument to attribute the high conductance peaks to a donor is provided by the valley splitting A. 
It is estimated from the dl jdVsd curves as the center of the cotunneling step, see data analysis in Subtraction of the 
Background section. The value of about 6.2 meV is indeed inconsistent with the typical valley splitting of < 1 meV 
of a Si/Si 02 quantum dolP^. 


Subtraction of the background 

In this section we show that the height of the Kondo peak is irrespective from the choice of the background model. 
The quantitative analysis of the Kondo contribution to conductance is performed by taking into account both inelastic 
cotunneling and thermal activated transport. By increasing microwave power the Kondo peak amplitude decreases, as 
a result of oscillating bias and decoherent microwave-induced processes like spin flip cotunneling; on the other hand, 
the background conductance in Coulomb blocked regions increases due to photon assisted tunnelling. The cotunneling 
step enhances as well, as predicted by Flensber^^ and experimentally confirmed in Refs. EH and I521 To discriminate 
the effects, the data analysis reported in Ref. ESI allows to isolate the Kondo physics from other competitive effects. 
In Figure the zero bias Kondo resonance for TV = 2 is shown at different microwave power values, namely the dark 
signal and at the nominal applied power -17 dBm, -15 dBm, -11 dBm and -8 dBm, before and after data analysis, with 
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■ Data 

N = 3 — Parabola 



Vsd (mV) 


FIG. 9. Examples of the background fitting using three different functions: a parabola (red), double exponential (blue) and 
double Boltzmann (green). The data refer to the N = 3 Coulomb diamond with a 30 GHz radiation except the top-left panel. 
At -15 dBm the Kondo resonance is completely suppressed. 


the microwave frequency fixed at 30 GHz. Without any background subtraction at the higher microwave power, i.e. 
from -11 to -8 dBm, it is difficult to distinguish the Kondo resonance from the Coulomb valley bottom. In Figurej^ 
the Kondo peak is decoupled from background according to the following procedure. We first systematically test three 
different background functions for = 1,2,3. A parabolic, a double exponential and a double Boltzmann fittings 
are imposed in the proximity of the resonance^^. As an example we report in Fig. |^the data analysis for N = 3. 
Generally parabola or a double exponential functions return a lower bound for experimental background signal, while 
the double Boltzmann gives a fit closer to data. Table [ll| reports the adjusted coefficients obtained by fitting the 
data of Figure It confirms the qualitative conclusions previously derived: the parabolic fit is characterized by the 
lowest among the tested functional forms, thus signaling the biggest deviation from the experimental data. The 
double exponential returns a R^ value close, but generally lower, to the double Boltzmann. Analogous results, not 
shown here, are obtained for = 1, 2. According to these results the data in Figure 5 of the main text and Figs. 


Parabola Double exponential Double Boltzmann 


0.982 

0.985 

0.985 

0.986 

0.987 

0.992 

0.980 

0.982 

0.982 

0.965 

0.981 

0.986 


TABLE II. : Adjusted for the fits shown in Fig. 


[^[^[^[^are extracted after a 2-Boltzmann background subtraction. Once determined the best functional form 
of the background signal, the next step is the data analysis of the resulting Kondo resonance. To evaluate the height 
of such resonances for N = 2,3 we use a simple Lorentzian (the theoretical shape of the Kondo peak at non-zero 
temperatures is still debated): 


2A w 


9 = go + 


TT 4(14d - 14 )^ + 


( 3 ) 
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Background: 2-Boltzmann 



Vsd (mV) 


FIG. 10. Example of a Lorentzian fit performed after background subtraction. The data shown are taken at a nominal 
microwave power of -19 dBm at 15 GHz, in the N = 2 valley. 


Background: 2-Boltzmann 



FIG. 11. Example of data analysis for = 1. The non-zero bias resonance is characterized by a non trivial shape, showing a 
shoulder at -6.2 mV due to the cotunneling step. 


where the fitting parameters are go, A, w and 14- The Lorentzian fits adequately the data, as reported in Fig. |10| 
where a -19 dBm irradiation is applied at 15 GHz. 

In order to evaluate and discriminate spin-valley Kondo effect from inelastic cotunneling at TV = 1 we adopt a modified 
approach. The cotunneling signal is predicted to be a step smoothed by temperature, superimposed on a background; 
in some cases it can be characterized by a cusp-like resonanc^^. In absence of an overall functional form describing 
the simultaneous presence of Kondo effect, inelastic cotunneling and background signal we choose to fit our data 
with the sum of two Lorentzians (one fits the Kondo resonance and one fits the cotunneling resonance) superimposed 
to each of the three fitting functions mentioned above. At 4.2 K and in dark conditions, the two Lorentzian result 
centered at Vsd = —6.2 mV and Vsd = —7 mV, which are identified as the cotunneling and Kondo peak respectively. 
The best E? (= 0.92) are obtained with the double Boltzmann choice, returning the satisfactory fits of Figure 
As shown in Fig. |I3| the Kondo resonance height behaviour is insensitive to the chosen background function. It is 
worthy to note that the fitting procedure with double Lorentzians returns good results when the two resonances are 
comparable. As soon as one prevails on the other this procedure will never return a complete suppression, producing 
an artificial saturation of the height of the decreasing peak. 

We finally remark with Figure [T^ that the Kondo peak decreases by increasing of microwave power whereas the 
cotunneling behaves oppositely. The net result consists in a shift of the position of the conductance maximum 
towards less negative biases, as visible in Figure [T^. 
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FIG. 12. (a) Contrary to the Kondo resonance, the cotunneling contribution to conductance is enhanced by lowering the 
temperature (b) Same experimental data as in Fig. shown after background subtraction at different powers of 30 GHz 
signal. The Kondo contribution, indeed, decreases whereas the cotunneling one increases resulting in a shift of the maximum 
of the peak. 



FIG. 13. Height of Kondo resonance for = 1 at 30 GHz as function of nominal power. 
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Background: 2-Boltzmann 



P (dBm) 


FIG. 14. Height of Kondo peaks for N = 2 as function of nominal microwave power at frequencies of 4, 15, 30 and 40 GHz. 
The power of the onset of suppression P* is conventionally chosen at the interception between the dark value and the linear 
decrease of the maxima. 


Rescaling procedure of nominal microwave powers 


As briefly reported in the main text, the transmission efficiency of the ac signal depends on several elements: the 
nominal 10 dB attenuator, the coaxial cable, the antenna, possible modes excited in the sample holder, behavior of 
the capping layers of the device. Each element of the line between the generator and the quantum system has its 
own frequency dependent impedance. We assume that at a fixed frequency the power of the microwave held delivered 
by the generator is proportional to the square of the amplitude of the signal seen by the quantum system: P oc V^. 
In order to determine the constant of proportionality, we follow the procedure devised by Elzerman et al.^^ once 
denoted as P* the nominal power of irradiation when the suppression of the Kondo peak begins, we set V* = ksTes/e, 
where Tea = 5.2 K from the thermal broadening of the first Coulomb peak at base temperature. Thus, since from our 
assumption P* jV*"^ = PjV^, we can link each nominal power P to the amplitude of the oscillations V^j applied to the 
sample. Our conventional onset of suppression P* is extracted by intercepting the dark value of the Kondo resonance 
with that of the trend when suppression is higher than about 15% (10 times the experimental uncertainty), see Fig. 
[T4| We have checked the robustness of such a procedure by subtracting the three different background functions, 
thus obtaining different V^. We have found no significant deviations, which demonstrates that the scaling behaviour 
observed is independent from our data analysis. Figure is an example for the = 1 case. The evaluation of 
the amplitude lA; from the nominal power P allows to replot the heights of the Kondo peaks as a function of the 
dimensionless parameter eV^/hv. In Figure [T^ we demonstrate that the procedure adopted to calibrate the frequency 
of the microwave signal is not dependent on the background subtraction. A comparison of the three panels of Fig. 
highlights some common features. With any background choice a collapse onto an unique curve is seen for the highest 
frequencies, i.e. 30 GHz and 40 GHz. At 4 GHz a weak dependence of the peak hight on the parameter eV^jhv 
is anyway present. The difference lies in the absolute value of this dimensionless parameter due to the different 
goodness of the fits. As discussed above we observe the best for the double Boltzmann, reasonable for the 
double exponential and quite acceptable results for the parabolic signal. The goodness of the fits influences the peak 
height evaluation, and by consequence the extraction of eV^jhiv. 


Microwave effect on a Kondo spin-valley shell 


In Figure we show an example of the results obtained perturbing the system with a microwave irradiation. Here 
the responses for A^ = 1,2,3 to a 30 GHz radiation are reported. At this frequency the microwave irradiation causes 
a suppression of the Kondo resonances for the three occupancies. 
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FIG. 15. Evolution of the Kondo resonance height at different microwave frequency in function of the parameter eV^jhv. In 
the left panel the subtracted background is fitted with a parabola, in the middle one it is used a double exponential, whereas 
in the right one a donble Boltzmann has been subtracted. 



FIG. 16. Kondo resonance suppression as a function of the nominal power P of the 30 GHz microwave field for the first 
spin-valley shell. 
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